Topological change and impedance spectrum of rat olfactory receptor 17: A 
comparative analysis with bovine rhodopsin and bacteriorhodopsin 
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We present a theoretical investigation on possible selection of olfactory receptors (ORs) as sensing 
components of nanobiosensors. Accordingly, we generate the impedance spectra of the rat OR 
17 in the native and activated state and analyze their differences. In this way, we connect the 
protein morphological transformation, caused by the sensing action, with its change of electrical 
impedance. The results are compared with those obtained by studying the best known protein of 
the GPCR family: bovine rhodopsin. Our investigations indicate that a change in morphology goes 
with a change in impedance spectrum mostly associated with a decrease of the static impedance 
up to about 60 % of the initial value, in qualitative agreement with existing experiments on rat 
OR 17. The predictiveness of the model is tested successfully for the case of recent experiments 
on bacteriorhodopsin. The present results point to a promising development of a new class of 
nanobiosensors based on the electrical properties of GPCR and other sensing proteins. 
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I. INTRODUCTION 

Nanobiosensors constitute a new frontier in the field of 
electronic devices, both for their reduced dimensions and 
for their high specialization in sensing action pj. One 
of the most promising challenge is producing more effi- 
cient devices owing to a direct use of biological sensing 
materials such as proteins. Sensing proteins are able to 
catch a specific molecule (ligand), e.g., odour, hormone, 
drug, light (photons), etc., then undergoing a change of 
their conformation. In vivo, the conformational change 
activates a chain of biochemical events that culminates 
in the transmission and detection of the captured infor- 
mation. In vitro, it is not easy to reproduce the entire 
chain and, probably, it is also not convenient. Thus, the 
capture of a specific ligand should be monitored looking 
only at some peculiar aspects of the capture mechanism, 
like the conformational change in the protein. 

In previous works 0, S] we described the possibil- 
ity to monitor a conformational change by using the 
change of the impedance spectrum associated with the 
single protein. Our model was based on the conjecture 
that if the protein is able to sustain a charge transport, 
then the macroscopic electrical response, specifically the 
impedance spectrum, would depend on the path followed 
by the charges and thus on the protein structure. When 
the structure changes, the impedance spectrum is ex- 
pected to change in a correlated way 0. 

In the present paper we apply the model to a cou- 
ple of sensing proteins, the 17 rat olfactory receptor Q 
and the bovine rhodopsin Both these proteins per- 
tain to the huge family of the seven-helices transmem- 
brane receptors, the so called G protein coupled receptors 
(GPCRs) p]. Therefore, they share a similar behavior in 
the sensing action: Both proteins undergo a conforma- 
tional change when going from the native to the activated 
state, then activating a G protein; the former when cap- 
turing a specific odorant, the latter when capturing a 



photon. As a consequence, it is potentially possible to 
reveal their conformational change by measuring some 
variations of their electrical properties. In this respect, it 
is crucial to determine whether this change is sufficiently 
large to be detectable within the experimental resolution, 
thus making possible to transduce the sensing action of 
the protein into an electrical signal. Recent experiments 
have related the sensing action of the rat OR 17 to the 
capture of a specific ligand. The OR 17 was deposited on 
a self- assembled multilayer (SAM) and the variation of its 
impedance spectrum was revealed by means of an electro- 
chemical impedance spectroscopy (EIS) characterization. 
These results disclose the possibility of using proteins as 
very refined sensors, able to determine, through an elec- 
trical signal, the presence and the concentration of the 
substance to be detected. 

The paper is organized as follows. Section II briefly 
recalls the theoretical model. Sections III and IV report 
the results together with their physical interpretation. 
Major conclusions are drawn in Sect. V. 



II. THEORY 

Our aim is to determine the variation of the protein 
impedance spectrum Q subsequent to a conformational 
change, by taking into account only geometrical changes 
in the tertiary structure of the protein. To this purpose, 
we assume that the 3D structures of the protein in its na- 
tive and activated state are known from the protein data 
bank (PDB) available in literature Q or similarly. In the 
present case, since the 3D structures are only partially 
known both for rhodopsin and OR 17, we made use of a 
homology modeling for both the proteins , as detailed 
later. Then we introduce a two step procedure. 

First, the protein in an assigned state is represented 
by means of a topological network (graph). Each node 
of this graph corresponds to a single amino acid: If two 
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amino acids are closer than an assigned cut-off distance, 
R , then a link is drawn between the corresponding nodes. 
The distance R definitely selects the connectivity of the 
network, and is here taken as an adjustable parameter 
of the model. Accordingly, the choice of its value has 
to be headed by some physical constraints and finally 
fixed by comparison with experiments. At this stage, 
we limit ourselves to observe that the network is built 
up with a strategy analog to that of random graphs, by 
using R , instead of a constant probability, to select the 
effective links [TTJ] . Therefore, for this model, we have 
the possibility to recover some of the results obtained in 
the context of random and complex networks (llj . 

Second, by associating an elemental impedance with 
each link, the graph becomes an impedance network. In 
other words, each link mimics a privileged channel for 
the charge transfer and/or charge polarization [H, ll2T |. 
By taking appropriate contacts and applying an exter- 
nal bias, the electrical network is solved within a linear 
circuit analysis. 

In a simple approach, the spatial position of each node 
is taken to correspond to the a-carbon of each amino 
acid. Each link connecting two nodes is replaced by an 
elementary impedance consisting of an Ohmic resistance 
connected in parallel with a parallel plate capacitor, filled 
with a dielectric. In such a way, we can describe both 
charge motion and electrical polarization. The elemen- 
tary impedance between nodes i and j, Z^j, is given ex- 
plicitly by [|[ : 
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where A% t j = n(R 2 — is the cross-sectional area 

between two spheres of radius R centered on the i-th and 
j-th node, respectively; Uj is the distance between these 
centers, p is the resistivity, taken to be the same for ev- 
ery amino acid, with the indicative value of p = 10 O 
m; i = \J — 1 is the imaginary unit, eo is the vacuum per- 
mittivity, and u> is the circular frequency of the applied 
voltage. The relative dielectric constant of the couple of 
i,j amino acids, ejj, is expressed in terms of the intrinsic 
polarizability of the i, j amino acids [l3[. 

By positioning the input and output electrical con- 
tacts, respectively, on the first and last node (correspond- 
ing to the first and last amino acid of the protein se- 
quential structure) , the network is solved within a linear 
Kirchhoff scheme and its global impedance spectrum is 
calculated within the convenient frequency range 0-100 
kHz. For a given interaction radius, the spectrum de- 
pends on the network structure. The change of the pro- 
tein state, by implying a change in the network structure, 
will lead to a different spectrum of the global impedance. 
Such a variation of the impedance spectrum is here taken 
as an estimate of the electrical response consequent the 
protein conformational change. Since the resistivity as- 
sociated with an amino acid is taken to be independent of 
geometry and applied voltage, the above estimate should 
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FIG. 1: Sketch of the backbone of bovine rhodopsin and rat 
OR 17 obtained with the homology modeling, a), b) refer to 
bovine rhodopsin in the native and activated states, respec- 
tively; c), d) refer to rat OR 17 in the native and activated 
states, respectively. 



be considered as a reference threshold value. By con- 
sidering more detailed microscopic mechanisms of charge 
14| . like tunneling, and/or the possible influ- 



transport 

ence of the medium surrounding the protein, the electri- 
cal sensitivity to the conformational change can be fur- 
ther modulated and, in particular, amplified. The role 
of the medium (lipids) surrounding the protein is here 
neglected. Indeed, the medium is not involved in the 
conformational change and its effect cancels when con- 
sidering the variation of the impedance as done here. 



A. Materials and methods 

The primary structure, i.e., the sequence of amino 
acids, is presently available for many proteins. However, 
the tertiary structure, i.e., the spatial organization, is 
known only for a very few of them and often is rather 
incomplete @. In particular, when the representations 
of different states of the protein are given, they mostly 
refer to different experimental conditions (X-ray resolu- 
tion, pH, temperature, etc). The body of all these fac- 
tors produces a large error bar in the assessment of the 
protein topology and in the possibility of faithfully dis- 
criminate different configurations Q. At present there 
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FIG. 2: Variation of the link total number vs the interaction 
radius for bovine rhodopsin and rat OR 17. The horizontal 
dashed line serves as a reference for negligible variations. 



are not complete representations of the tertiary struc- 
ture for both the native and activated state of bovine 
rhodopsin, taken in the same experimental conditions. 
Even worse is the state of art of rat OR 17, for which 
only the primary sequence is taken by experiments. At 
present, for both the proteins, only the tertiary structure 
of the native states has been predicted, with first prin- 
ciple methods (l5| . For a reliable modeling, we require 
that both the native and activated 3D structures origi- 
nate by the same experimental or theoretical conditions. 
For this purpose, the tertiary structures of native and ac- 
tivated states for rhodopsin and OR 17 were obtained by 
means of the MODELLER software The template 
PDB entries were the same for both the proteins: 1JFP 
and 1LN6 for the native and activated state, respectively, 
complemented with 1F88, chain A. Both 1JFP and 1LN6 
were produced in the same experimental conditions. 

Figures 1 (a)-(d) report the backbone of bovine 
rhodopsin and rat OR 17 in their natural and activated 
states, as determined by the homology modeling. 

The method we adopt for comparing the change in 
impedance spectrum associated with the change of con- 
figuration is the analysis in the corresponding Nyquist 
plot. This plot is obtained by drawing the negative imag- 
inary part versus the real part of the global impedance, 
within a given frequency range (typ ically from 1 mHz 
to 100 kHz as in experiments For a single RC 

parallel circuit, the shape of the Nyquist plot appears as 
a perfect semi-circle. In the present case, owing to the 
presence of many different RC circuits, the plots slightly 
depart from the perfect semi-circle shape, as more gen- 
erally described in terms of Cole-Cole plots [3. [l7j. 



III. RESULTS AND DISCUSSION 

The topological network obtained from the protein 
mapping is in general dependent on the value of R . In- 
deed, by increasing R , the degree of each node, i.e., the 
number of links spanned by each node fill ] , increases in 
a way that reflects the node distribution in the network. 
For very large values of R , the node degree saturates to 
the expected value N — 1, with N the number of nodes. 
Consequently, the total link number shows a character- 
istic sigmoid-likc profile, which saturates to the value 
N(N- l)/2. 

Figure 2 reports the relative variation of the total link 
number associated with the native and the activated 
state as function of R , for rat OR 17 and rhodopsin, re- 
spectively. For small values of R , e.g., less than 6 A , the 
network is poorly connected, which implies a large sen- 
sitivity of the number of links on R . On the other hand, 
for large values of R , e.g., larger than about 20 A , most 
of the distances between nodes are smaller than R and 
quite all the nodes are connected: It implies a smooth 
variation of the total number of links. At these asymp- 
totic values of R , differences in the relative distance be- 
tween amino acids, induced by a conformational change, 
cannot be appreciated by the model. In the intermedi- 
ate region of R values, the number of links strictly de- 
pends on the distance between amino acids and thus the 
model makes possible to resolve the different states of the 
protein. Accordingly, the values of R useful for making 
the network able to capture the variation of impedance 
subsequent to a conformational change, should be cho- 
sen in the intermediate region 6-20 A , which also over- 
laps with the range of physical interactions among amino 

acids dGlE!. 

Figure 3 reports the degree distribution for relevant 
values of the interaction radius. Here we observe that, 
for small values of R , say 6 A , the distribution is rather 
peaked implying that each amino acid interacts practi- 
cally with the same number of nearest neighbors, inde- 
pendently of its position in the protein. The distribution 
broads at increasing values of R, say 9-12 A, where, in 
any case, several small but detectable peaks appear in 
the distribution profile. These peaks signal the differ- 
ent clustering of the protein main structures (helices and 
sheets). Accordingly, for the values of R in the range 
6-12 A , the network is maximally able to monitor the 
internal structures of the protein. When R is further in- 
creased, e.g., R larger than 12 A , the degree distribution 
becomes very broadened and no longer able to provide 
information on the internal structures of the protein. 
Finally, in the case the conformational change mainly in- 
volves displacements of nearest neighbors, the topological 
network has to be constructed by choosing values of R in 
the range 6-9 A . This choice leads to quite different pat- 
terns for the natural and activated structure. Otherwise, 
when the conformational change involves displacements 
of entire structures like helices, the values of the inter- 
action radius should be chosen in a slightly wider range, 
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FIG. 3: Degree distribution for the native states of bovine 
rhodopsin and rat OR 17 (in the insert). The curves corre- 
spond to values of the interaction radius in the range from 6 
to 12 A . Distributions obtained with larger values of R loose 
of resolution in this scale. 




FIG. 4: Relative variation of the static impedance of bovine 
rhodopsin and rat OR 17 vs the interaction radius 



say 9-12 A . 

Figure 4 reports the relative variation of the static 
impedance (resistance) as function of the interacting ra- 
dius for rhodopsin and rat OR 17. We can observe that 
the maximal variation is in the range 5-20 A , for both the 
proteins. This is in accordance with what shown in Figs. 
2 and 3. Furthermore, from this figure it emerges that the 
variation of resistance is never smaller than 10 % for both 
the proteins and, as a consequence, their conformational 
change is large enough to be revealed from variations of 
their electrical properties, in particular, from EIS mea- 
surements. Wc notice that the variation of resistance 
is significantly larger than the link difference, indicat- 
ing the greater relevance of the specificity (length and 
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FIG. 5: Nyquist plots for the rhodopsin [frames (a) and (c)] 
and rat OR 17 [frames (b) and (d)]. Both the real and imag- 
inary parts of impedance have been normalized to the value 
of the static impedance of the corresponding activated state 
of bovine rhodopsin and rat OR 17, respectively. Continu- 
ous lines refer to the native state, and dashed lines to the 
activated state, respectively. 



position) of the links rather than of their number. This 
outcome downsizes the importance of the precise value of 
R : In fact, the electrical resolution of the conformational 
change is practicable within a wide range of R values. 
Therefore, the qualitative accordance with experimental 
data is preserved to a large extent, while the quantitative 
agreement can be achieved by fine tuning the value of R . 

On this basis, to evidence the most relevant features 
that a conformational change in the protein can induce on 
its impedance spectrum, we take two R values, 6 A and 12 
A for both rhodopsin and rat OR 17. Figure 5 compares 
the Nyquist plots for the native and activated states of 
rhodopsin [Figs. 5 (a) and (c) in the left column] and of 
rat OR 17 [Figs. 5 (b) and (d) in the right column]. For 
both proteins the shape of the plots is close to a semi- 
circle, apart a small squeezing for R = 6 A . This result 
indicates that a smaller distribution of RC circuits is re- 
flected in a deviation from the otherwise perfect semicir- 
cle shape, as expected Q . An increasing resolution with 
increasing values of R is evidenced for both proteins, with 
the best value obtained for R = 12 A : This result con- 
firms the presence of a large conformational change 0, G3 
In particular, for rat OR 17 the resistance between the 
configuration is found to exhibit a significative larger dif- 
ference than for rhodopsin at both the considered values 
of R . Finally, we can notice that while in rhodopsin the 
native state shows a smaller value of the static impedance 
with respect to the activated state, in OR 17 it happens 
the opposite. The origin of this difference can be traced 
back to the shrinking trend of the activated structure of 
OR 17, and in the opposite behavior of rhodopsin, as can 
be guessed from Fig. 1. The theoretical results for rat 
OR 17 compare well with recent EIS experiments carried 
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out on OR 17 samples deposited on SAM in the pres- 
ence of specific odorants, heptanal and octanal Q- The 
experimental Nyquist spectra exhibited a near semicircle 
shape in the region of intermediate low frequencies with a 
decrease of the resistance for about 20 % in the presence 
of octanal at a concentration of 10~ 4 mol/1, in qualita- 
tive agreement with the expectation of the present model. 
The experiments gave evidence also of a Warburg resis- 
tance at the smallest frequency, here neglected. At this 
stage, since we are comparing macroscopic measurements 
with a microscopic model (single protein) the qualitative 
agreement found is considered to be satisfactory by con- 
firming, in essence, that the capture of an odorant by an 
OR belonging to the GPCR family can be monitored by 
means of electrical measurements. 




(R-RJ(A) 



A. Global properties 

We conclude this section by reporting the behavior of 
global protein conductance and capacitance, as functions 
of R . Figure 6 shows the static conductance, G = l/Z(0), 
versus the reduced cut-off radius ( R - R o ) , both for the na- 
tive and the activated state of rhodopsin, where R less 
than Ro «3.6 A is the value of R, which corresponds 
to zero conductance (i.e., disconnected network). The 
static conductance is found to follow a power law behav- 
ior: G cx (R — Ro) 7 , where the exponent 7 takes two dif- 
ferent values for two different regions of R values: 7 ~ 2 
for(R-Ro) >30A,and 7 ~ 3for(R-R ) < 20 A (see 
Fig. 6), for both the native and activated states. We 
can explain this result by observing that from Eq. ([1]) 
it emerges a R 2 -dependence of the single conductance. 
Then, for increasing number of links, the parallel connec- 
tions increase more than those in series and this produces 
a further R -dependence. When the number of links ap- 
proaches the saturation value, only the R 2 -dependence 
survives. As expected, the global capacitance is found to 
follow the same behavior of the conductance, as shown 
in Fig. 7. Analogous results have been obtained, but not 
shown, for rat OR 17. 



IV. A FURTHER TEST OF THE NETWORK 
MODEL 

As a further validation of the network model, we have 
considered the case of bacteriorhodopsin. The change of 
its electrical properties due to the conformational change 
induced by green light has been recently investigated 
through current voltage measurements pol . [2l| . Bac- 
teriorhodopsin, i.e., a seven helices light receptor like 
rhodopsin, is present in archaea instead of mammals 
and works with a different mechanism of signal trans- 
mission. Like rhodopsin, its light-sensitive part is the 
retinal, which is located deep inside the protein. When 
retinal captures photons, it changes its shape and in- 
duces the protein conformational change. However, in 



FIG. 6: Static conductance of the bovine rhodopsin pro- 
tein. The symbols refer to calculations; the straight lines are 
the best fit curves. The continuous line refers to the native 
rhodopsin (Rho) and the dashed line refer to the activated 
rhodopsin (Meta II). 




FIG. 7: Global capacitance of the bovine rhodopsin protein. 
The symbols refer to calculation; the straight lines are the best 
fit curves. The continuous line refers to the native rhodopsin 
(Rho) and the dashed line refer to the activated rhodopsin 
(Meta II). 



bacteriorhodopsin retinal changes its shape from straight 
to bent, just the opposite of what happens in rhodopsin 
@ . Accordingly, for bacteriorhodopsin we expect a pat- 
tern of the Nyquist plot associated with the transfor- 
mation, opposite to that exhibited by bovine rhodopsin, 
see Fig. 5, where the activated state exhibits a static 
impedance greater than that of the native state. To 
this purpose, Fig. 8 reports the Nyquist plot for two 
different couples of native and activated representations 
of bacteriorhodopsin, 1FBB (native) - 1FBK (activated) 
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FIG. 8: Nyquist plots for two couples of native and acti- 
vated states of bacteriorhodopsin: 1FBB - 1FBK and 2NTU 
- 2NTW, native and activated states of the protein, respec- 
tively. Plots have been obtained with R =6 A . 

and 2NTU (native) - 2NTW (activated), for R = 6 A. 
The plots confirm that the static impedance of the na- 
tive state is significantly greater (up to 15 %) than that 
of the activated state, which is, in fact, a behavior op- 
posite to that of bovine rhodopsin. On the other hand, 
this result conforms with the expectations of the topolog- 
ical analysis and agrees with the experimental evidence 
of an increase of conductance when passing from the na- 
tive (i.e,. in the dark) to the activated state (i.e., in the 
presence of green light) [2(| ■ These findings are taken as 
a further validation of the network model here reported. 



their function. In fact, the model allows an estimation 
of the magnitude of electrical impedance variation asso- 
ciated with the activation of the protein. These kind of 
measurements, at present quite difficult to be carried out, 
are, anyway, of great interest. In fact, they can disclose 
the possibility to produce a new generation of sensors of 
nanometric dimensions and sensible to very small quanti- 
ties of hormones, neurotransmitters, toxic products, and 
so on 

When applied to rhodopsin and rat OR 17, the 
impedance network predicts the possibility to resolve the 
different configurations pertaining to the native and ac- 
tivated state, by the difference exhibited by the corre- 
sponding Nyquist plots. The resolution level is predicted 
to achieve a maximum of about 60 % in the variation 
of the static impedance. EIS of rat OR 17 sensing ac- 
tion has evidenced changes of the impedance spectrum 
up to about 20 % when specific odorants have been 
introduced on the samples anchoring the ORs. These re- 
sults are taken as an indirect validation of the present 
modeling. As a further validation test, the methodology 
here proposed was also applied to the bacteriorhodopsin. 
The model correctly predicts the opposite behavior of 
the static impedance of this protein compared to that 
of bovine rhodopsin. Also this behavior is in qualitative 
agreement with experimental measurements [20j. As fi- 
nal remark, we mention the possible role played by the 
surrounding media of the protein in determining the ab- 
solute value of the impedance. Here we neglected this 
effect, which we expect to be rather insensitive to the 
conformational change. 
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V. CONCLUSIONS 

We have investigate the impedance spectrum of rat 
OR 17 and carried out a comparative analysis with that 
of the homologous GPCR, bovine rhodopsin. The used 
model is a powerful tool in the analysis of the confor- 
mational changes these proteins undergo when exerting 
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